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Abstract 

Lab-on-chip devices and point-of-care diagnostic chip devices are composed of many different 
components such as nanosensors that must be able to communicate with other components within 
the device. Molecular communication is a promising solution for on-chip communication. In particular, 
kinesin driven microtubule (MT) motility is an effective means of transferring information particles from 
one component to another. However, finding an optimal shape for these channels can be challenging. 

In this paper we derive a mathematical optimization model that can be used to find the optimal channel 
shape and dimensions for any transmission period. We derive three specific models for the rectangular 
channels, regular polygonal channels, and regular polygonal ring channels. We show that the optimal 
channel shapes are the square-shaped channel for the rectangular channel, and circular-shaped channel 
for the other classes of shapes. Finally, we show that among all 2 dimensional shapes the optimal design 
choice that maximizes information rate is the circular-shaped channel. 
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I. Introduction 

With the advancements in the field of nanotechnology, applications such as lab-on-chip devices 
and point-of-care diagnostic chips Q are becoming a reality. In most of these applications 
communication between different components in each device, or between a biological entity and 
a component in the device, such as a nanosensor, is required. Therefore, one of the obstacles 
that must be overcome before many of these applications can be fully realized, is devising a 
communication system for small scales Q, 0. 

Inspired by nature, one of the most promising solutions is molecular communication 0, 0, 
where molecules released by a transmitter device propagate in a fluidic environment and transfer 
information to a receiver device. Information can be conveyed by encoding messages into the 
release timing 0, number 0, 0, concentration p^ , or identities of particles pT] |. Different 
types of propagation are possible for transporting information particles in on-chip molecular 


communication such as: diffusion [12|-[14|, diffusion with flow [15|, active transport using 
molecular motors and cytoskeletal filaments p^ , bacterial assisted propagation pT] |, pT| , p0 , 
and kinesin molecular motors moving over immobilized microtubule (MT) tracks p9| , p0| . 
Molecular communication is biocompatible, and can be very energy efficient. Because of these 


properties, molecular communication has attracted a lot of attention in recent years |211. It has 


many potential applications in biomedical engineering from lab-on-chip devices p2| and point 
of-care diagnostic devices 


to targeted drug delivery p4| . However, most previous works 
have considered diffusion based molecular communication, which can be very slow for on-chip 
applications. Moreover, in 0 it was shown that active transport using MT filaments gliding 
over kinesin covered substrate can generate higher information rates compared to flow based 
propagation over larger separations distances. Other types of propagations are either difficult to 


implement or have not been implemented yet. Moreover, in [25|-[27| it is shown that electrical 
currents can be used to control the speed and direction of the MTs. This makes kinesin-driven MT 
motility a suitable form of propagation for on-chip molecular communication and transport, with 
potential applications in lab-on-chip devices, point-of-care diagnostic devices, organ-on-a-chip 
devices, and microfluid devices. 

For the kinesin-microtubule based molecular communication to be fully realized, guide posts 
for design and implementation are required. For example, in microfluidics it is well known that 
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different ehannel designs eould be used to implement different funetionalities in droplet mieroflu- 


idies [281. However, finding the optimal design for these ehannels ean be diffieult beeause of the 
system eomplexities sueh as the random motion of the MTs. Moreover, experimental trial-and- 
error-based design is very laborious and time-eonsuming. Therefore, simulation environments 


based on experimentally validated models sueh as [29|, [30| ean be a very effeetive tool for 
providing design guide posts. For example in our previous works Q, [ [3T] |, optimization of 
the shape of the transmission zone (i.e. loeation plaeement of transmitters) is eonsidered using 
eomputer simulations, and it was shown that optimal transmission zone is along of the walls of 


the ehannel. This is beeause usually MTs follow the ehannel walls [32|. Moreover, in [33|, |34| 
we have shown that the shape of the ehannel ean have a huge effeet on information transmission 
rate, and we provide basie guidepost for designing the shape of the ehannel. 

In this work, we extend the preliminary results presented in our previous works and derive 
an optimization model that eould be used to find the optimal ehannel shape and dimensions. 
Our results ean be used in two ways. First, in on-ehip transport applieations, our results eould 
be used to maximize the transport rate. On-ehip transport is a very important proeess in lab- 
on-ehip deviees. For example, in [ |35| , kinesin-mierotubule transport seheme is used for high- 
throughput moleeular transport and assembly. Our results ean also be used in on-ehip moleeular 
eommunieation applieations to maximize the ehannel eapaeity of the system. We show that using 
the optimized ehannel shape ean inerease the eapaeity by more than 1.5 bits per ehannel use and 
in some eases by more than 2.5 bits per ehannel use. This ean be a signiheant gain in sequential 
data transmission. 

To setup our optimization formula, first we derive a eomplete model relating the number of 
MT trips during a single transmission period and the shape of the ehannel. This is our first major 
eontribution. We then use this model to find the optimal shape and dimensions for reetangular 
ehannels, regular polygonal ehannels, and regular polygonal ring-shaped ehannels. This is our 
seeond major eontribution. Finally, we show that for all possible two dimensional ehannel shapes, 
the optimal ehannel is the eireular-shaped ehannel. Moreover, for a partieular transmission rate 
we show that we eould use our optimization model to find the optimal radius for this eirele. 
This is our third major eontribution. We verify all our results using eomputer simulations. 

The rest of this paper is organized as follows. In Seetion an overview of the kinesin- 
mierotubule based moleeular eommunieation system is presented. A general information trans- 
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Fig. 1: Depiction of the communication environment. 


mission model is presented in Section]^ In Section |IV] we derive the channel shape optimization 
model, and in Section |V] we use the derived model to find the optimal channel design strategies. 


The results are then compared with computer simulations in Section VI Concluding remarks 


are given in Section VII 


II. Kinesin-Microtubule based Molecuear Communication Channels 

For our model we assume that an area in the kinesin covered channel is designated as the 
transmission zone, where the transmitter and the information particles released by the transmitter 
reside. Similarly, an area in the channel is designated as the receiver zone, where the receiver and 
its receptors are placed. Generally, we assume the transmission zone and the receiver zone are 
located at the opposite ends of the channel. Moreover, the information particles released by the 
transmitter remain stationary at the transmission zone until they are picked up by MT filaments 
gliding over kinesin covered substrate. The MT filaments then carry these information particles 
to the receiver zone, where they are unloaded and delivered to the receiver. The information can 
be encoded in the type, number, release timing, and concentration of the particles released and 
delivered. Moreover, the information particles can be encapsulated inside liposomes to protect 
them from denaturation (e.g., molecular deformation and cleavage caused by enzymatic attacks 
or changes in pH of the outer aqueous phase) in the propagation environment. 

For the loading and unloading information particles, we assume deoxyribonucleic acid (DNA) 


hybridization bonds are employed as shown in [36|. In this scheme, information particles, MT 
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Fig. 2: Different shape elasses and their parameters. 


filaments, transmission zone, and the reeeiver zones are eovered with single stranded DNAs 
(ssDNA). When information particles are released by the transmitter, they anchor themselves 
to the transmission zone through DNA hybridization bonds until an MT filament passes at the 
vicinity of the anchored particle. At this point the particle is loaded on to the passing MT again 
through DNA hybridization. When a loaded MT enters this receiver zone, the loaded particles 
are unloaded using yet another DNA hybridization bond. This process is summarized in Figure 
[^and the reader is referred to [361 for detailed explanation and laboratory demonstration of this 
technique. 


A. Shape of the Channel 

In this paper we assume the channel can be formed to have a wide variety of different shapes, 
and we find the optimal shape that maximizes information transmission rate among all possible 
shapes. In particular, we assume that channel shape can belong to any of the following three 
classes of shapes: rectangular, regular polygon, and regular polygon ring. Regular polygons, are 
equiangular (all angles are equal in measure) and equilateral (all sides have the same length), 
and they include a large class of geometric shapes, ranging from equilateral triangles, squares, 
pentagons, and hexagons, all the way up to a circle as number of sides approaches infinity. 
Figure shows an example from each of the three different shape classes. In practice, these 
channels, as well as channel having any other shape can be created using similar procedures 
used in [ |7^ , [ |30| |, [ [3^ . 


Depending on the class of the channel shape, different parameters can be used to further 











6 


adjust the shape within the shape class. For rectangular channels two parameters can be used to 
adjust the dimensions and the shape: width w and length 1. When both the width and the length 
are equal in value the channel would be a square channel. Similarly, for regular polygons two 
parameters can be used to further adjust the shape: the number of sides n, and the radius of the 
circumscribed circle r (this radius can also be defined as the distance from the center point of the 
regular polygon to one of the edges of the polygon). For ring-shaped channels three parameters 
are used to define the shape: the number of sides n, the radius of the inner circumscribed circle 
Tj, and the radius of the outer circumscribed circle Tq. Figure summarizes all these parameters 
for different classes of shapes. 

Regardless of the shape of the channel, without loss of generality, we assume the transmission 
zone is always on the left side of the channel and the receiver zone is along the right side of the 
channel. In [[^, it was shown that the optimal transmission zone is along the walls of the channel 
since MTs mostly glide very close to channel walls. This positioning increases the chance of 
an MT picking up an information particle during its trips. Therefore, in this work we always 
assume the transmission zone is along the walls of the channel regardless of the channel shape. 
Moreover, the minimum separation distance between the start of the transmission zone and the 
receiver zone are always the same in our numerical evaluations of Section |V^ regardless of the 
shape of the channel. 


B. Simulating the Channel 


In [161, [361 the authors show that creating kinesin-microtubule based channels described in 
the previous section is possible in wet labs. However, studying the resulting molecular communi¬ 
cation channel, and solving design and optimization problems using laboratory experimentation 
is extremely difficult because these experiment are time consuming, laborious, and expensive. 
To overcome this issue, computer simulations have been used in previous works [|30|. 


Similarly, in this work we use Monte-Carlo techniques proposed in [29| to simulate the motion 
of the MTs over a kinesin covered substrate. The equations for the motion of the MT, which 


were developed and verified experimentally in [29|, are given below. 


The motion of the MT is largely regular, although the effects of Brownian motion cause 
random fluctuations. Since the MTs move only in the x-y directions, and do not move in the z 
direction (along the height of the channel), we consider a two-dimensional simulation of MTs 
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for At time intervals. Given some initial position (xo,yo) at time t = 0, for any integer k > 0, 
the motion of the MT is given by the sequenee of eoordinates (xi, yi) for i = 1,2,.Eaeh 
eoordinate {xi,yi) represents the position of the MT’s head at the end of the time t = iAt, 
where 


Xi = Xi-i + Ar cos 9i, 

(1) 

yi = yi-i + Ar sin 9i. 

(2) 


In this ease, the step size Ar at eaeh step is an independent and identieally distributed (i.i.d.) 
Gaussian random variable with mean and varianee 


E[Ar] = UavgAf, 

(3) 

Var[Ar] = 2DAt, 

(4) 


where Uavg is the average veloeity of the MT, and D is the MT’s diffusion eoeffieient. The angle 
9i is no longer independent from step to step: instead, for some step-to-step angular ehange A9, 
we have that 

9i = A9 + 9i_,. (5) 


Now, for eaeh step, A9 is an i.i.d. Gaussian-distributed random variable with mean and varianee 

E[A9] = 0, (6) 

Var[A«] = (7) 


where Lp is the persistenee length of the MT’s trajeetory. Following [29|, in ease of a eollision 
with a boundary, we assume that the MT sets 9i so as to follow the boundary. Finally, our 
simulator ean model any polygon-shaped ehannel. 

We also use the grid loading meehanism proposed in Q, to simulate the information partieles’ 
loading and unloading. For loading an information partiele, the MT filament must drive elose to 
the anehored partiele. Therefore, we divide the transmission zone into a square grid, where the 
length of eaeh square in the grid is the same as the diameter of the partieles. We then distribute 
partieles randomly and uniformly among the squares in the grid. In general, we assume that 
the MTs ean load multiple partieles, whieh we know to be possible based on lab experiments; 
thus, if an MT enters a square whieh is oeeupied by a partiele, and it has an empty loading 





Fig. 3: An example of the simulation environment. 


slot available, we assume the MT loads that partiele. For unloading, we assume all the loaded 
partieles are unloaded as soon as an MT enters the reeeiver zone. 

Figure shows a sample simulation. In this figure, the ehannel is a regular polygon ring with 
parameters n = 8, = 20 |am, and Tq = 25 |am. Other parameters used for simulation are: 

simulation time steps of AT = 0.1 seeonds, MT diffusion eoeffieient D = 2.0 ■ 10“^ p.m^/s, 
average speed of the MT Vavg = 0.5 p.m/s, and persistenee length of the MT trajeetory Lp = 111 
p.m. We also assume the size of the information partieles is 1 |am, the average length of the MTs 
is 10p.m, and eaeh MT ean load up to 5 information partieles in one trip from the transmission 
zone to the reeeiver zone. These parameters are all seleeted based on experimental observations 
of DNA eovered MTs moving over a kinesin eovered substrate, and the reader is referred to [j^ 
for more details about the simulation environment. In the rest of the paper we will always use 
these parameters unless it is speeified otherwise. 

In Figure the MT initially starts on top of the ehannel. It will then move towards the 
transmission zone (blue/eyan grid) along the left wall of the ehannel. The MT loads the first 5 
information partieles on its way and eontinues to the reeeiver zone where it unloads the loaded 
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information particles. Notice that although the MT enters grid squares that contain information 
particles, after 5 loaded particles, it does not load anything else because the maximum load 
capacity is reached. 


III. Information Transmission Model 

In molecular communication, messages can be encoded into information carrying particles 
using many different schemes. Information can be encoded into the number, concentration, or 
type of the particles released. One key observation in molecular communication is that, regardless 
of the encoding scheme, information is always transmitted through mass transfer (i.e. transfer 
of particles). Therefore, in this section we present a general information transmission model for 
mass transfer, which is independent of the information encoding scheme. 

Let X = - jTmax} be the set of possible particles that could be released by the 

transmitter, where Xmax is the maximum number of particles the transmitter can release. Let 
X G X he the number of information particles released into the medium by the transmitter. 
We define time per channel use (TPCU) as a predefined amount of time T representing the 
duration of a single transmission session. Because of the random motion of the MT filaments, 
given this time limit, information might not be perfectly conveyed to the receiver. For example, 
it is possible that some of the particles will not arrive at the receiver after T has elapsed and 
therefore there is some information loss. This effect is similar in nature to the noise introduced 
by conventional electronic or electromagnetic channels. 

Let G X represent the number of information particles that arrive at the destination 
after TPCU duration T. Because the receiver does not have any prior knowledge of the number 
of particles released by the transmitter, at the receiver X is a discrete random variable given 
by probability mass function (PMF) fx{x). Similarly, given X particles were released by the 
transmitter, and the TPCU duration T, is a discrete random variable given by conditional 
PMF I x). This conditional PMF is very important because it characterizes the 

channel completely. It can be used to calculate important channel parameters such as channel 


capacity [371, the maximum rate at which any communication system can reliably transmit 
information over a noisy channel. 

Assuming that the transmitter and receiver are perfect (i.e. they release and receive the 
information particles perfectly), and that the information particles are not lost in the environment. 
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the only factor that effects the PMF \ x) is the random motion of the MTs. Although 

this assumption may not hold in practice, by assuming perfect transmitter and receiver, we only 
focus on the effects of channel-shape on the information rate. Moreover, for some applications 
such as pure mass transport in lab-on-chip devices, we are only concerned with mass transport as 
opposed to a communication system. The PMF /y(T)|x(|/^^^ | x) could be used as an important 
performance measure in both cases. For a given TPCU duration T, as the number of MT trips 
between the transmitter and the receiver increases, the number of particles that could potentially 
be transported increases. This results in higher mass transfer and potentially higher achievable 
information rates. Based on this observation we develop an optimization model which can be 
used to find the optimal channel shape. 

IV. Channel Shape Optimization Model 

Let the set Q be the set of all possible cross sectional geometric shapes the channel could 
have, where the cross sectional shape is the shape of the channel when viewed from the top. 
In this paper we fist assume this set contains all the rectangular, regular polygonal, and regular 
polygonal ring shapes and then extend the result to all two dimensional shapes. Our goal is to 
find the optimal shape in these sets, given system parameters such as TPCU, concentration of 
MTs, and average speed of the MTs. To setup this optimization problem, we model the effects 
of channel shape on the number of MT trips. 

As discussed in the previous section, the number of MT trips during the TPCU interval has a 
direct effect on transportation of information particles. Without loss of generality let G ^ be 
the channel under consideration. Therefore, let K^'^\g) be the number of trips during the TPCU 
interval T inside the channel shape g, where a single MT trip is defined as the movement of 
the MT from anywhere in the channel to the transmission zone and then the receiver zone. For 
example, a single MT trip is shown in Figure]^ After the MT completes its first trip, subsequent 
trips are defined as the movement of the MT from the receiver zone to the transmission zone and 
back. During any trip, an MT can deliver 0 or more information particles (up to its maximum 
load capacity). Because the motion of the MT filaments are random in nature, K^'^'>{g) is a 
random variable. Therefore, to derive our optimization model we use the expected value of 

Let the random variable K^\g) be the number of trips for a single MT during the TPCU 
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Fig. 4: Average number of trips approximation eompared with Monte Carlo simulations. 

interval, T. Let Uavg be the average speed of the MTs, P{g) be the perimeter of the ehannel 
shape g G Q. A good estimate for the average number of MT trips, when a single MT is inside 
the channel is given by 



( 8 ) 


This approximation is based on the observation that on average from each trip, the MT travels 
a distance equivalent to the perimeter of the channel. This assumption is verified using our 
Monte Carlo simulation environment. The results of the simulations are shown in Figure In 
this figure, we have considered a sample shape from each of our three different shape classes, 
and have shown that the approximations in Equation Q (solid lines) are a good estimate of the 
actual average number of MT trips (point plots). 

In practice, there are typically more than one MT inside the channel. Furthermore, the number 
of MTs is not a constant, but is dependent on the volume of fluid in the channel: the DNA- 
covered MTs are prepared in chemical solutions, and therefore have a constant concentration 
inside the solution. Let C be the concentration of MTs as number of MTs per unit volume, 
and h be the height of the channel, and ^(^f) be the cross sectional area of the channel with 
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geometric shape g E Q. Then the number of MTs inside the channel is given by 

M{g) = A{g) X hx C, 


( 9 ) 


where M{g) is the number of MTs inside channel g. 

For each MT in the channel let the random variables Kf"\g) be the number of trips during 
T seconds inside the channel with shape g, for i = 1, 2, • • • , M. Then, the total number of MT 
trips during T seconds is given by 


M 




( 10 ) 


i=l 


The average number of MT trips during T seconds is therefore calculated as 


M 




( 11 ) 


2=1 


We assume that the number of trips for individual MTs are independent and identically dis¬ 
tributed, because they are chemically similar and don’t interact with each other. Generally, this 


assumption has been found to be valid in laboratory experiments p^ , p0| . Due to the identical 
distributions, the equation simplifies to 


M 


S|A''"’'(9)1 = E A|AT’(9)1 = M(g)E[KP{g)]. 


( 12 ) 


2=1 


Using the approximation shown in Equation Q, and Equation Q, the total number of trips can 
be estimated as 


E[K^^\g)] ^Txv^^gXCxhx 


PigV 


(13) 


where A{g) and P{g) are the cross sectional area and perimeter of the channel with geometric 
shape g E Q, respectively. The channel shape optimization problem can then be formulated as 


uiSiX E[K^P {g)]. 
g&Q 


(14) 


Assuming all the other parameters are constant except the cross sectional shape (including the 
height of the channel), the optimization problem becomes 

'm' 


max 


g&g [P{g)_ 


(15) 


Equation (15) is very interesting because it states that the optimal shape is the one with the 
largest area to perimeter ratio. However, there is an important constraint that must be satisfied. 
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The perimeter of the ehannel must be small enough such that a single MT can complete at least 
a single trip. Therefore, we assume the perimeter must be such that its length can be traveled 
by the MT, on average, during the TPCU interval T. This constraint ensures that the perimeter 
is small enough such that on average enough information particles can be delivered during the 
given TPCU duration. Finally, the optimization problem can be written in its complete form as 



subject to P{g) < Twavg- 


(16) 


V. Optimal Shape Analysis 


Using the optimization formula derived in the previous section, we analyze the optimal channel 
design for each shape class individually, and then the optimal design over all. 

1) Rectangular Channels: Let ^Rec be the set of all rectangular shapes. Because any rectan¬ 
gular shape can be characterized by the two parameters width w and length I, the Equation ( [T^ 
becomes 


w X I 


max - - 

{'W^l)GQ^ec 2x0 “t” 2/ 


subject to 2w + 21 < Twavg- 


(17) 


Solving this optimization problem, the optimal channel design given both TPCU interval T and 
average speed of the MTs Vavg, is w = I = 0.25Tx;avg. Therefore, for rectangular channels, the 
optimal channel shape is always the square shaped channel. 

2) Regular Polygonal Channels: Let ^poiy be the set of all regular polygons. Regular polygons 
can be characterized by two parameters: the number of sides n, and the radius of the circum¬ 
scribed circle r. Note that the set of all square shapes are also contained within the set ^poiy 
(i.e. when n = 4). Using these parameters Equation ([T^ becomes 


O.Snr^ sin(27r/n) 


max 

(n,r)e5poiy 


2nr sin(7r/?7,) 


subject to 2nrsin(7r/n) < Twavg. 


(18) 
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Simplifying this equation we get 


max 0.5r cos(7r/n) 
(n,r) eCpoiy 


subjeet to nr sin(7r/?7,) < 



(19) 


where we have used the faet that sin(2n) = 2sin(n) cos(m). Based on this equation the optimal 
ehannel is the one with n = oo and r = Travg/2vr. Therefore, when regular polygons are 
eonsidered the optimal shape is the eireular-shaped ehannels. Beeause squares are also in the set 
of all regular polygons, we eonelude that eireular-shaped ehannels are also better than square¬ 
shaped ehannels. 

3) Regular Polygonal Ring Channels: Let ^Ring be the set of all polygonal ring shapes. The 
elements of this set ean be eharaeterized using three parameters: the number sides n, the radius of 
the inner polygon’s eireumseribed eirele r*, and the radius of the outer polygon’s eireumseribed 
eirele Tq. Note that ^Ring D ^poiy, whieh means that the set ^Ring eontains all the regular polygonal 
shapes (this is aehieved by setting r* = 0). Using these parameters the optimization problem 
beeomes 


max 

in,ri,ro)&QRiug 


0.5n{r‘^ — rf) sin(27r/?7,) 
2nro sin(7r/n) 


subjeet to 2nroSin(7r/n) < Tuavg- 


( 20 ) 


Simplifying this equation we get 


max 

(n,ri,ro)e6Rmg 


0.5 


^2 _ ^2 


subjeet to nroSm{Tr/n) < 


cos(7r/n) 
Tu. 


avg 


( 21 ) 


Solving Equation 21. the optimal ehannel has parameters n = oo. To = Tuavg/2vr, and r* = 0. 


This means the optimal ehannel is the eireular-shaped ehannel. This may seem surprising at 
first, sinee one would expeet that a ring-shaped eireular ehannel would be better than a solid 
eireular-shaped ehannel at transporting information partieles. However, beeause typieally MTs 
follow the walls of the ehannel, and that the number of MTs in a ehannel is proportional to the 


volume of the ehannel [32|, it beeomes apparent that the solid eireular ehannel would be better 


than a ring-shaped ehannel. 
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4) Overall Optimal Channel Design: From the solution to the optimization formulas for eaeh 
of the three different shape elasses, it is apparent that the optimal ehannel shape is the eireular- 
shaped ehannel when the set of geometrie shapes ^rpr = ^Rec U ^Poiy U ^Ring is eonsidered. 
Moreover, based on Equation ( [T^ the optimal ehannel shape must have the largest area to 
perimeter ratio. Therefore, if the perimeter of the ehannel is fixed sueh that the eonstraint in the 
optimization is satisfied with an equality, the optimal ehannel shape would be a eirele for all 
two dimensional geometrie shapes. 


VI. Results and Discussions 

To verify our optimization formulas and their results, we rely on Monte Carlo simulation. For 


these simulations we use the same parameters diseussed in Seetion II-B Moreover, we assume 
the height of the ehannel is always a eonstant = 10 p.m regardless of the eross seetional shape 
of the ehannel. The eoneentration of the MTs is also assumed to be (7 = 0.001 MT/fL unless 
speeified otherwise. The number of MT in the ehannel is always assumed to be 


MT=[A{g)hC\, 


( 22 ) 


unless speeified otherwise. 


For the performanee measure we rely on channel capacity [37|, whieh is the maximum rate 
at whieh information ean be transmitted reliably. Channel eapaeity is given by 


C = max /(X; 

fx{x) 


(23) 


where is the mutual information between X and Mutual information is defined 


as 


/(X; V('^)) = E 


logs 




(24) 


where fymixiy^'^^lx) represents the probability of reeeiving symbol at the destination 
during TPCU interval T, given that symbol x was transmitted by the souree; fx{x) represents 
the probability of transmitting symbol x; and E[-] represents expeetation. Using Monte Carlo 
simulations the PMF /y(t)|x(?/*'^^|t) is estimated for a given ehannel shape. For more details 
about the Monte Carlo simulations for estimating this PMF please see Q, and for analytieal 
estimation of this PMF using Markov ehain models please refer to [j^. In this work, we use 
Monte Carlo simulations sinee they are more aeeurate in multi-MT environments. Using the 
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TPCU = 160s, x^^ = 40 


TPCU = 240 s, = 40 



(a) 


(b) 


Fig. 5: Channel capacities of different rectangular channel shapes. The channel with the highest 
capacity is shown using the blue dot. For TPCU of 160 seconds the optimal channel is 20 p.m 
X 20 |am (a), and for TPCU value of 240 seconds the optimal channel dimensions is 30 p.m x 
30 p.m (b). 


estimated PMF /y(T)|x(|/*'^^|ai), the maximizing PMF fxix) and hence the channel capacity is 


calculated using Blahut-Arimoto algorithm [39|, 

We now consider each of three shape classes discussed in the previous sections, and use the 
Equations ( [T7] ), ( [T9| ), and ( |2T] ) to find the optimal channel dimensions for each case. These results 
are then compared to the channel capacities obtained from Monte Carlo simulations. 


A. Rectangular Channels 

First, we consider rectangular channels. We consider two values of TPCU: T = 160 seconds, 
and T = 240 seconds. Using the solution to Equation ( [T7| ), the optimal channels are square 
channels with dimensions w = I = 0.25 x 160 x 0.5 = 20 qm and w = I = 0.25 x 240 x 0.5 = 30 
qm for the T = 160 and T = 240, respectively. To verify this result, we simulate all the 
rectangular channels with widths and lengths raging from 20 qm to 50 qm in 5 qm steps. 
We assume Xmax (he- the maximum number of particles the transmitter can release) is 40. The 
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TPCU value of 320 seconds 



Fig. 6: Channel capacity of regular polygonal channels with constant perimeter of 160 p.m 
versus maximum number of particles that can be released at the transmitter. 


channel capacity is then calculated based on the simulations. Figure shows the results for 
T = 160 seconds (a), and T = 240 seconds (b). As can be seen the channels with the largest 
channel capacity are the 20 p.m x 20 p.m, and 30 p.m x 30 p.m, respectively. These results are 
in perfect agreement with the results obtained from our optimization model. 

B. Regular Polygonal Channels 

In Section |V] we showed that for regular polygonal channels the optimal shape is a circular¬ 
shaped channel. To verify this result, we first consider a constant TPCU value of 320 seconds. 
From Equation ( [T^ , the maximum channel perimeter that satisfies the constraint is 160 p.m. 
Therefore, we assume that the channel perimeter is fixed at this value and show that the capacity 
increases a the channel shape become more circular by increasing the number of sides n. In 
particular we consider the following channels: square channel of length 40 qm (equivalent to 
a polygon with n = 4), and regular polygonal channels with parameters (n = 8, r = 26.13 
qm), (n = 12, r = 25.76 qm), (n = 16, r = 25.63 qm), and (n = 20, r = 25.57 qm). 
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TPCU value of 250 seconds 



Fig. 7: Channel capacity of four circular-shaped channels versus maximum number of particles 
that can be released at the transmitter. 


Based on our optimization formula we expect the 20-sided channel to be the optimal among all 
these channels since it is more circular. Figure shows the channel capacity for each of these 
channels obtained through Monte Carlo simulations versus the maximum number of particles the 
transmitter can release Xmax- As can be seen, the optimal channel is indeed the 20-sided regular 
polygon. Moreover, from the obtained pattern it is obvious that as the number of sides increase 
the channel capacity increases. This supports the results obtained from our optimization formula 
that circular channels are optimal. 


To evaluate the performance of Equation (19) we consider the TPCU interval of T = 250 
seconds. The solution of this optimization problem is a channel with parameters n = oo and 
r Ri 20 qm. Because in our simulation environment the number of sides must be a finite number, 
we use the value of n = 20 to simulate a perfect circle. We then simulate four different regular 
polygonal channels with parameters: (n = 20, r = 17 p.m), (n = 20, r = 20 qm), {n = 20, 
r = 22.75 qm), and (n = 20, r = 25.57 qm). Figure shows the channel capacity of each 
channel versus the maximum number of particles the transmitter could release. As can be seen 
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TPCU value of 320 seconds 



Fig. 8: Channel capacity of ring-shaped channels versus maximum number of particles that can 
be released at the transmitter. 


in the figure, the channel with the radius of 20 p.m achieves the highest capacity among the four 
channels. This result is in agreement with the optimal solution derived from our optimization 
formula. 


C. Regular Polygonal Ring Channels 

In this section, we validate our optimization formula for polygonal ring-shaped channels 


presented in Equation (21) through simulations. We assume the TPCU interval is fixed at 320 
seconds. Using the solution to Equation ( |2T] ), the optimal channel dimensions are n = oo, 
Ti = 0, and To 25.57 qm. In Eigurej^we have already shown that the solid circular-shaped 
channel is optimal among regular polygons. Therefore, we only consider ring-shaped channels 
for validation. In particular, we consider the following ring-shaped channels: (n = 20, To = 25.57 
qm, Tj = 0 qm), (n = 20, Vg = 25.57 qm, = 10 qm), {n = 20, Vg = 25.57 qm, r* = 15 
qm), and (n = 20, Vg = 25.57 qm, r* = 20 qm). Eigure shows the channel capacity of these 
channels versus the maximum number of particles the transmitter could release. We can see 
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from the figure that a solid eireular ehannel aehieves the highest eapaeity as predieted by our 
optimization formula. 

Based on all these results, we ean see that a solid eireular-shaped ehannel is the optimal shape 
that maximizes information rate, when stationary kinesin with mobile MT filaments are used for 
transportation in moleeular eommunieation systems. Moreover, as the number of sides inerease 
and the ehannel beeomes eloser to eireular, the performanee gains of adding more sides deereases. 
Based on the obtained results, there is negligible performanee differenee between ehannels with 
the number of sides greater than 20. Therefore, these ehannels are suffieient representation of 
perfeetly eireular-shaped ehannels. 


VII. Conclusion 

In this work we eonsidered the problem of finding the optimal ehannel shape for on-ehip 
moleeular eommunieation based on kinesin driven mierotubule (MT) motility. To setup our 
optimization model, we first derived a mathematieal model relating the average number of MT 
trips to the shape of the ehannel. Using this model we then presented the three optimization 
formulas and their respeetive solutions for three different elasses of shapes: the reetangular, 
regular polygonal, and regular polygonal rings. We showed that the optimal solution is the 
square ehannel for the reetangular ehannel shapes, and eireular shape for the regular polygonal 
and ring-shaped ehannels. Finally, we showed that when all elasses of 2 dimensional shaped are 
eonsidered the eireular-shaped ehannel tend to be the one that aehieves the highest information 
rate. Moreover, using our optimization formula we showed that the dimensions of the optimal 
ehannel ean be reliably estimated. 
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